The PDF of fluid particle acceleration in turbulent flow with underlying normal 

distribution of velocity fluctuations 
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We describe a formal procedure to obtain and specify the general form of a marginal distribution 
for the Lagrangian acceleration of fluid particle in developed turbulent flow using Langevin type 
equation and the assumption that velocity fluctuation u follows a normal distribution with zero 
mean, in accord to the Heisenberg-Yaglom picture. For a particular representation, (3 = exp[«], of 
the fluctuating parameter /3, we reproduce the underlying log-normal distribution and the associated 
marginal distribution, which was found to be in a very good agreement with the new experimental 
data by Crawford, Mordant, and Bodenschatz on the acceleration statistics. We discuss on arising 
possibilities to make refinements of the log-normal model. 
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I. INTRODUCTION 



In the context of Tsallis formalism Q and general- 
ized statistics approach @. [J3 to developed turbulence, 
C. Beck [E| used the underlying log-normal distribution, 
/(/?), to describe fluctuation of the parameter (3 entering 
some Langevin equation for the acceleration a of fluid 
particle in the Lagrangian frame; see also 0,0- For the 
case of a linear drift force, the resulting marginal proba- 
bility density function of the acceleration is [5| 
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exp 
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(1) 



where m = exp[s 2 ] provides a unit variance and s is a fit- 
ting parameter, s 2 = 3.0. This distribution was found 
to be in a very good agreement with the recent La- 
grangian experimental data by Porta, Voth, Crawford, 
Alexander, and Bodenschatz p|, the new data by Craw- 
ford, Mordant, and Bodenschatz (the Taylor microscale 
Reynolds number is R\ — 690, the normalized acceler- 
ation range is [—60, 60] 3 a, the Kolmogorov timescale 
t v is resolved) [jg, Mordant, Delour, Leveque, Arneodo, 
and Pinton (R\ = 740, a <G [-20, 20], t v is not resolved) 
|10|. and direct numerical simulations of the Navier- 
Stokes equations by Kraichnan and Gotoh (R\ = 380, 
a e [-150,150]) Ej. 

The same approach with the underlying % 2 distribu- 
tion, /(/?), of fluctuations resulted in an analytically ex- 
plicit distribution fl^. 



P(a) 
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where C is a normalization constant, Vq = 4, and q = 3/2 
(Tsallis entropic index) are due to the theory, and a c is 
a free parameter used for a fitting, a c — 36.0. This dis- 
tribution was found to be in a good agreement with the 
experiments Q|. The result is based on the consider- 
ation of Ref . [Lj ■ 

In general, the marginal distribution is defined as 



P(a) 



df3 P(a\0)f(J3), 



(3) 



where P(a\(3) is a probability density function associated 
to a surrogate dynamical equation, the Langevin equa- 
tion for the acceleration a |l3( , 



d t a = ~/F(a) + aL(t). 



(4) 



Here, F(a) is a "drift force", L(t) is delta-correlated 
Langevin source (Gaussian white noise, (L(t)) = 0, 
{L(t)L(t')) = 2S(t-t')), 13 = -f/a 2 is assumed to be a fluc- 
tuating real positive parameter. For constant parameters 
7 and a, this model assumes that the stochastic process 
Q is Markovian, and P(a\/3) is found as a stationary 
solution of the associated Fokker-Planck equation. For a 
linear drift force, F(a) — —a, the stationary conditional 
distribution, P(a\(3), is found to be 



P(a|/3) = C(/3)exp[-/3a 2 /2], 



(5) 



where C(j3) is a normalization constant, a € [—00,00]. 

The function /(/?) entering Eq. (3) is probability den- 
sity function arising from the assumption that the pa- 
rameter f3 is a stochastic variable. For a constant /?, 
the Gaussian solution JSJl meets the large scale (large 
time increment) statistics but it fails to describe observed 
Reynolds number dependent stretched exponential tails 
of the experimental small scale P(a) [£| that correspond 
to high probability to find extremely big values of the 
fluid particle acceleration in the developed turbulent flow. 

The interest in studying Langevin type equations to 
describe developed turbulence is motivated by the re- 
cent high precision Lagrangian experiments ja 0, Il0| - 
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which give an important dynamical information and new 
look to the intermittency in fluid turbulence. Time re- 
sponse characteristics of the tracer polystyrene particle 
and the precision in the experiments |a, |9j allow to re- 
solve about 1/20 of the Kolmogorov time and 1/40 of the 
Kolmogorov length (R\ — 970) so that the acceleration 
is really resolved and particle follows rare violent events 
(within 7% of the ideal value of acceleration even at the 
high Reynolds numbers studied there), with the collected 
statistical data sufficient to establish convergence of the 
fourth moment of acceleration. 

In Sec. 2, we construct a simple model that allows one 
to derive the marginal probability density function, P(a), 
of the acceleration of fluid particle in turbulent flow, with 
underlying normally distributed velocity fluctuations. In 
Sec. 3, we reproduce the result of Ref. Q as a particu- 
lar case, and propose a generalization of the log-normal 
model. In Sec. 4, we briefly summarize the results and 
make a few remarks. 



II. THE MARGINAL DISTRIBUTION WITH 
UNDERLYING NORMALLY DISTRIBUTED 
VELOCITY FLUCTUATIONS 

While the conditional distribution P(a\f3) can be found 
starting from the Langevin equation (gj, we do not have 
sufficiently strong theoretical requirements to determine 
a unique form of the distribution f{(3) except for that it 
should obey general conditions, such as that the integral 
in Eq. should converge and the resulting marginal dis- 
tribution P(a) should be normalizable. To select a par- 
ticular form of the distribution /(/3), one can use ad hoc 
statistical distributions or make assumptions stemming 
from considerations of the developed turbulence. 

As a first step, one can associate the parameter f3 with 
the mean energy dissipation rate e, and use the relation- 
ship of e with the fluid velocity Vi , within the framework 
of Kolmogorov scaling theory of developed turbulence. 

Second, as shown by Renner, Peinke, and Fried- 
rich 01 > the small scale intermittency can be traced back 
to a stochastic nature of the averaged energy dissipation 
rate 0] , and the Markovian condition for the stochastic 
process is fulfilled at scales larger than the Taylor mi- 
croscale A. Thus, a stochastic nature of [5 can be related 
to the stochastic energy dissipation rate (measured in 
Eulerian experiments), which we denote by e. 

For example, the assumptions (3 ~ e and e ~ u 1 > dis- 
cussed at length by Beck Jjj| (we do not repeat it here for 
brevity), for the fluctuating, averaged over a ball of radius 
r, energy dissipation e and velocity fluctuations Ui at the 
Kolmogorov scale were used to propose (3 ~ ^ uf and to 
select the % 2 distribution of j3. Three-dimensionality of 
the space was used explicitly to determine the value of 
the free parameter (Tsallis entropic index) in this model. 
A comparison of the formalism based on Eq. J3J with the 
Sawford model 16] implies (3 = 2C ~ 2 a i/ 1 / 2 e~ 3 / 2 @ so 
that with the replacement of e by the fluctuating e and 



using again e ~ uf , one can obtain a different relation, 
P ~ X>r 3 - The Kolmo gorov's 1941 approach implies 
e ~ it 3 , where u stands for the rms velocity, which can 
also be used to try to relate e to the velocity fluctuations. 
Although successful in capturing main features of the ex- 
perimental data, different models use different powers of 
e to represent j3 that makes a theoretical problem when 
selecting an appropriate distribution of (3. 

It is however common to the above representations that 
velocity fluctuations, Ui, are assumed to be normally dis- 
tributed with zero mean, and that statistics of (3, e, and 
Ui are related to each other due to some functional de- 
pendencies between their characteristics holding in the 
inertial range. A direct dependence of statistical prop- 
erties of the acceleration a on velocity fluctuations was 
established in the wellknown Heisenberg-Yaglom theory. 
This gives a different look to the problem since the La- 
grangian description refers to individual trajectories of 
fluid particles as compared with the well established Eu- 
lerian framework, in which the intermittency is under- 
stood in terms of anomalous scaling of the moments of ve- 
locity increments in space related to the stochastic nature 
of the energy dissipation rate (non-dynamical descrip- 
tion). From a simplified dynamical point of view, one 
can think of that the tracer particle trajectory extends 
to a large region (few integral length scales, in the exper- 
iments) crossing during the course subregions character- 
ized by different local amplitudes of velocity fluctuations 
which are randomly distributed in space and varies with 
time (Brownian like motion of the particle driven by the 
stochastic delta-correlated force with stochastically slow 
varying intensity). The Lagrangian velocity autocorre- 
lation function is known to decay very slowly, to vanish 
at times bigger than the integral time scale. This view 
is valid for sufficiently high Reynolds numbers and meets 
that provided by the stochastic energy dissipation rate 
argued earlier [13j . This is also in an agreement with the 
remark made in Ref. Q that extremely high accelerations 
seem to be associated with coherent structures (geomet- 
rically these may be thought of as very thin tornadoes) 
which persist for many Kolmogorov times, substantially 
longer than the correlation time of the acceleration com- 
ponents. At lower Reynolds numbers, R\ < 500, the 
acceleration is increasingly coupled to large scales of the 
flow in a dynamical way so that a different approach may 
be required to describe such a situation to a good accu- 
racy. 

In the Lagrangian frame, it is then natural to relate pa- 
rameters of the stochastic dynamical equation (0} , which 
describes the acceleration in the context of a general- 
ized Brownian motion, to velocity fluctuations due to the 
Heisenberg-Yaglom picture of developed turbulence. Ul- 
timately of course we deal with a statistical description 
of the acceleration and the velocity fluctuations within 
the framework of Fokker-Planck equation associated to 
Eq. Q. Guided by this observation, one can develop a 
general consideration along this line of reasoning. 

In the present paper, we formulate some general re- 
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quirements for /(/?) using the assumption that the ve- 
locity fluctuations itj are Gaussian distributed random 
variables with zero means, (itj) = 0, and variances Sj, 



g(u t ) = C(s t )exp[-^/(2 S 2 )] 



(6) 



where C(sj) is normalization constant; Uj G [—00,00], 
i = 1,2,3. This assumption meets the Lagrangian ex- 
perimental data 8], for each component of the velocity 
fluctuations. 

It should be stressed that our treatment is made re- 
gardless to a particular functional dependence of (3 on e, 
or a particular form of /(/?), and provides a simple way to 
formalize the general model in accord to the Heisenberg- 
Yaglom picture, which relates some statistical properties 
of the acceleration to that of the velocity fluctuations. 

We assume that (3 depends on the velocity Vi the dy- 
namics of which is decoupled from that of the acceleration 
so that the well known fluctuating character of Vi for the 
turbulent flow determines fluctuations of (3. Below we 
drop the index i to simplify notation. Also, we treat the 
developed turbulent flow to be statistically isotropic and 
translation invariant thus discarding skewness effects; the 
statistical anisotropy is known to become very small with 
the increase of the Reynolds number, e.g., at R\ > 900. 

In the mean field approximation, we can represent v = 
(v) + u, where the fluctuation u is characterized by zero 
mean, (u) = 0. Hence, we can write 



(3 = p c ((v))+B((v),u) ) 



(7) 



where we formally separated out f3 c which depends only 
on (v). 

We suppose that the function B(u) = B((v),u) satis- 
fies the following natural requirements: (i) It is a suffi- 
ciently smooth invertible function with respect to u so 
that the inverse function u((3 — (3 C ) can be found; (ii) it 
maps the interval [—00,00] to [— /3 c ,oo], to provide pos- 
itiveness of (3; also one can restrict consideration by the 
requirement that (iii) there is one-to-one correspondence 
between (3 and u. 

In general, for (3 to be a stochastic variable the function 
B(u) should be Borel function of the stochastic variable 
u. We assume that u has an absolutely continuous dis- 
tribution function and the associated probability density 
function, in particular, that given by Eq. JSJ). Not any 
Borel function B(u) provides existence of a probability 
density function for (3. However, a probability density 
function for (3 a priori exists if B(u) is a monotonic func- 
tion. 

Using the general relation 



g(u)du = g{u(f3- (3 C )) 



du 



dft 



(8) 



and Eq. ©, we can make the identification f((3) 
g(u(f3-(3 c ))\du/df3\,i.e., 



/(/?) = C( S )exp 



[u(J3 - 0c 



2s 2 





du 




d(3 



(9) 



where C(s) is a normalization constant. The inverse 
function u{(3 — (3 C ) should provide positive definitcness 
and normalizability of the above function f{(3). From 
Eq. we then finally have 



P(a) = C(s) d(3P(a\(3)exp 



W-(3 c )f 



2s 2 





du 




dp 



(10) 



This equation allows one to calculate the distribution 
P(a) for given stationary solution P(a\(3) associated to 
the Langevin equation Q and the invertible (monotonic) 
function B{u). By construction, the structure of this 
model is such that the velocity fluctuation u, underlying 
the dynamics, is normally distributed with zero mean 
and variance s. Note that only the amplitude of u con- 
tributes P(a). Eq. ljrrj|> can be viewed as a specific case 
of the general equation 10 which captures one of the well 
known features of the turbulent flow, and thus enables to 
rule out some ad hoc types of distribution of (3. 

Note that in Eq. (|10f) we did not change the vari- 
able over which the integral is performed as compared 
to Eq. (PI). The main idea was to introduce the vari- 
able u which follows Gaussian distribution and encode 
various underlying distributions in the functional depen- 
dence (3 = B(u). This significantly reduces the class 
of admissible probability density functions f{(3) that we 
consider as one of the main results of our approach. 

In the next Section, we study a specific choice of B(u) 
which is relevant both from the mathematical and phys- 
ical points of view. 



III. THE LOG-NORMAL MODEL AND 
BEYOND 

As mentioned by Beck selection of the log- normal 
distribution f{J3) provides that for any power law, (3 ~ 
e K , the function ln/3 = relne guarantees that the func- 
tional form of the log-normal distribution does not 
change, which is viewed as a hint for a physical rele- 
vance of a log-normally distributed (3. Also, it is known 
that the probability density function of the averaged en- 
ergy dissipation rate is log-normal in agreement with Kol- 
mogorov's assumption in K62 theory. 

The exponential dependence, 



B{u) 



W u 



(11) 



where uq is constant, and (3 C = 0, provides a relevant 
example for which we get 



u = Uq ln/3, 



du 
W 



(12) 



The distribution @ then becomes the log-normal distri- 
bution, and Eq. (fTTTf) for P(a\(3) given by Eq. © repro- 
duces the marginal distribution Q), with uq = 1. 
Also, supposing (3 ~ e K , we obtain the relation 



a u/K 



(13) 
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where we put uq = 1 for simplicity. Again, one can see 
that the only measurable effect of the parameter k is that 
it scales the variance, s — > s/n, in Eq. lfT71|l. 

One can make the following refinement of the dynam- 
ical equation (0} using the relation (fill) , which corre- 
sponds to the log-normal model |QJ. The defining feature 
of the considered model is the assumption of slow fluc- 
tuating character of the composite parameter j3 = 7/c 2 - 
It can be easily shown that a noisy character of the drift 
coefficient 7, with F(a) = —a and delta-correlated Gaus- 
sian white noise 7 with non-zero mean, can be used to 
describe intermittency effects (power law tails) because it 
is related to wcllknown random multiplicative processes 
extensively studied in the literature. Below, we focus on 
a stochastic nature of the additive noise amplitude a and 
study its contribution to the intermittency separately. A 
joint effect of both the noisy 7 and random intensity of 
the additive noise a is out of scope of the present paper 
and can be studied elsewhere. 

From Eq. JTTJ} it follows a 2 = e r u / u ° so that the 
Langevin equation (0J becomes 

d t a = jF{a) + exp [-u/{2u )\ L(t), (14) 

where the velocity fluctuation u is treated statisti- 
cally and follows Gaussian distribution with zero mean, 
F(a) = —a, and 7 is taken to be constant to simplify 
consideration of the contribution of the additive noise in- 
tensity. The Langevin model i|14|) meets the Lagrangian 
experiments to a very good accuracy as it implies the 
probability density function of the form . 

Equation (|14|) assumes that the dynamics of La- 
grangian acceleration a (the acceleration in the comov- 
ing frame, = dfVi + v^dkVi), which comes mainly from 
small scales, and that of the fluctuating energy dissipa- 
tion rate e, or the related amplitude of velocity fluctua- 
tions, Ui = Vi — (vi), which is associated mainly to large 
scales, are weakly related to each other at high Reynolds 
numbers, R\ > 500, so that in the lowest approxima- 
tion it is natural to take them independent. At small 
time scale, the velocity fluctuation is characterized by 
approximately constant amplitude while its directional 
part changes wildly (cf. |8j). However, in the statistical 
context, which reflects Lagrangian dynamics at all the 
time scales, the acceleration is related to the amplitude 
of velocity fluctuation as it follows, for example, from the 
Heisenberg-Yaglom scaling, (a 2 ) ~ atiu?l 2 v~ x l 2 L~^l 2 , 
where u is the rms velocity and L is the integral scale 
length. This longstanding universal u 9 ^ 2 scaling was con- 
firmed by the recent Lagrangian experiments Q to a very 
high accuracy, for R\ > 500. The large scale origin of 
the additive noise stochastic intensity, slowly varying in 
time, is precisely the reason of the stretching of exponen- 
tial tails of the acceleration probability density function 
P(a). The observed Reynolds number dependence of the 
stretching can be thus qualitatively understood as the re- 
sult of increasing coupling of the acceleration dynamics 
to large scale dynamics at lower Reynolds numbers. 

In summary, the dynamical equation i|14fl reflects the 



above point of view by incorporating exponential of the 
dynamically independent u, with the prescribed exter- 
nal statistics, as the intensity of the additive noise. A 
full dynamical treatment of Eq. I|14l) with u viewed as a 
random function of time can be made elsewhere. Here, 
we note only that Langevin type equation with the ad- 
ditive noise of the form e"'*- 1 L(t), with w(t) a longtime 
correlated stochastic process, was recently proposed 
to describe the Lagrangian intermittency in the context 
of new multifractal random walk model by Bacry and 
Muzy, a continuous extension of discrete cascade models. 
Our approach (|14l) meets this result. The difference is 
that we approximate the random process u) to be triv- 
ially stationary and it obeys normal statistics. 

In accord to Eqs. (J5J and 114(1. the associated station- 
ary distribution for F(a) = —a is given by 

P(a\u) = C{u) exp[-e tl /"°a 2 /2], (15) 

where C(u) is a normalization constant. From this point 
of view, small deviations of the theoretical P( a) g iven by 
Eq. from the experimental data [j| (see |l2fl for de- 
tails) could be attributed to small deviations from Gaus- 
sian distribution of the velocity fluctuation u; the flatness 
was measured to be approximately 2.8 for the axial and 
3.2 for the longitudinal velocity component |8j compared 
with the flatness value 3 for an exact Gaussian distribu- 
tion. Eq. I|15fl can be used to account for such deviations 
in a phcnomcnological way, 

duC{u)exp[-e u/uo a 2 /2] gi (u), (16) 

-00 

where gi(u) is the distribution which is approximately 
Gaussian one. However, it should be noted that small 
departure from the experimental data can also arise from 
shortcomings of the considered one-dimensional Langevin 
model of developed turbulence. 

While it is evident that the three-dimensional Navier- 
Stokes equation with a Gaussian random forcing belongs 
to a class of non-linear stochastic equations for which 
one can associate some generalized Fokker-Planck equa- 
tions, it is a theoretical challenge to make a link between 
the Navier-Stokes equation and such a surrogate one- 
dimensional Langevin model, which is, of course, far from 
being a full model of the Lagrangian dynamics of fluid 
in the turbulent regime. This problem is addressed in 
the forthcoming paper [T?|. Here, we note only that the 
Lagrangian description simplifies the problem of finding 
the stationary probability measure since the advection 
term, which dominates in the developed turbulent flow, 
is incorporated to the definition of Lagrangian accelera- 
tion. Also, large scale dynamics and an interaction be- 
tween the large scale and small scale dynamics within 
the inertial range are important to provide understand- 
ing of the origin of noises (multiplicative and additive 
ones) used in such Langevin models and a weak depen- 
dence of the noise statistics on the small scale dynamics 
at high Reynolds numbers. 
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IV. CONCLUSIONS 

(i) The presented approach implies a specific form of 
the probability density function f((3), given by Eq. JjJJl, 
which stems from the assumption that velocity fluctua- 
tions are normally distributed with zero mean, Eq. JSJ, 
as it is confirmed by the experimental data 0,0 for high- 
Reynolds-number turbulent flows, to a good accuracy. 

(ii) Given the exponential dependence ((Til the pre- 
sented formalism leads to a consideration of the log- 
normal distribution /(/3), which was proved to be rel- 
evant from the experimental point of view ,5] . Also, such 
a distribution is in agreement with the log-normal distri- 
bution of mean energy dissipation in Kolmogorov 1962 
picture. An exponential relation of the type l|llfl re- 
quires a physical treatment in the context of Langevin, 
or Fokker-Planck, equation which can be made elsewhere. 
Here, we note that the exponential dependence on the ve- 
locity fluctuation, exp[w], indicates a specific variance of 
the absolute value of velocity increment during the time- 
evolution, and can be thought of as a kind of Lyapunov 



instability in the velocity space. 

(iii) One may also be interested in using some other 
distributions, instead of the normal one ©, to derive 
f((3) and the associated theoretical distribution P(a). 

(iv) Finally, one can try other appropriate functions 
B(u), for instance, 

B( M )= e "/ tl0 ]T Cfe u fe , (17) 

where Ck are constants, instead of to reproduce the 
experimental data on the acceleration statistics of fluid 
particle in turbulent flow. In particular, for (3 C = and 
the function B(u) = u 2 the map doubly covers the in- 
terval [0, oo] 9 (3, with each covering being a monotonic 
map, and we obtain u = ±^ffi and \du/d(3\ = P' 1 / 2 . The 
basic equation @ yields f(/3) ~ /?~ 1/2 exp[-/?/(2s 2 )], 
which is x 2 probability density function of order one. 
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